A curvature-driven effective attraction in multicomponent membranes 
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We study closed liquid membranes that segregate into three phases due to differences in the chemical and 
physical properties of its components. The shape and in-plane membrane arrangement of the phases are coupled 
through phase-specific bending energies and line tensions. We use simulated annealing Monte Carlo simulations 
to find low-energy structures, allowing both phase arrangement and membrane shape to relax. The three-phase 
system is the simplest one in which there are multiple interface pairs, allowing us to analyze interfacial prefer- 
ences and pairwise distinct line tensions. We observe the system's preference for interface pairs that maximize 
differences in spontaneous curvature. From a pattern selection perspective, this acts as an effective attraction 
between phases of most disparate spontaneous curvature. We show that this effective attraction is robust enough 
to persist even when the interface between these phases is the most penalized by line tension. This effect is 
driven by geometry and not by any explicit component-component interaction. 
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I. INTRODUCTION 

Multicomponent liquid membranes are pervasive in na- 
ture, and their phase behavior is important in many biolog- 
ical processes. Such membranes are of interest not only as 
constituents of complex biological systems but also as rich 
pattern-forming systems in themselves. They provide a fertile 
ground for exploring the role of geometry in surface pattern 
formation at the nanoscale, which is an important aspect of 
functional materials. They have been studied as basic mod- 
els for biological systems |0J and as promising candidates 
in the rational design of biocompatible materials [3]. Their 
promise is due in large part to their versatile phase behav- 
ior. In a liquid membrane, the membrane's constituents are 
free to rearrange. Although mixing is entropically favorable, 
this freedom has energetic consequences when there are mul- 
tiple kinds of membrane components, since different com- 
ponents may have different chemical and physical properties 
[J]. When the system adapts its arrangement to its conditions, 
varying shapes and component patterns can result. 

Patterns of segregated domains have been directly observed 
in experimentally created multicomponent giant unilamel- 
lar vesicles Id UiOII . These experiments agree with theoret- 
ical results in which patterns arise in response to diffe ring 
bending properties of the membrane's segregated phases [JTT|— 
[lrjll . These bending properties can include bending rigidities, 
saddle-splay moduli, and spontaneous (preferred) curvatures. 
When these properties are phase-specific, the phase arrange- 
ment is coupled to the membrane shape. 

A key aspect of surface pattern selection is the system's re- 
sponse to geometric constraint. Constraints include fixed sur- 
face area, since the energy required to compress a liquid mem- 
brane is orders of magnitude larger than the energy to bend it, 
and closure, since tears or holes in the membrane would lead 
to energetically unfavorable exposure of hydrophobic tails to 
surrounding water. For a closed membrane that does not ex- 
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change material with its environment, there may also be a 
constraint on enclosed volume. This does not apply in all 
cases; for example, lipid bilayers have been observed to adjust 
their volume by forming transient pores [17]. We use a con- 
stant volume assumption here for the sake of simplicity and 
note that relaxing this constraint does not qualitatively affect 
reported results. Under these constraints, a multicomponent 
membrane may not be able to fully indulge the preferences of 
all of its constituents simultaneously. 

Thus a strongly-segregated membrane may face this frus- 
tration: on the one hand, immiscibility (realized as a line 
tension between differing component domains) favors macro- 
scopic segregation and interfaces with minimal length ifTill ; on 
the other hand, when subject to geometric constraints, com- 
ponent bending preferences may be more favorably accom- 
modated by arrangements with extended interfaces, including 
arrangements with multiple domains. In this sense it is the 
matter of the length and location of interfaces that becomes 
the battlefront between these competing preferences. 

Most studies so far have focused on two-phase systems. In 
two-phase systems there is only one phase pair, hence only 
one type of interface. Introducing a third phase is a nontrivial 
extension. The presence of a third phase provides three phase 
pairs and thus introduces the crucial feature of interfacial pref- 
erences: the three-phase system has freedom not only over 
length and location of interfaces but also over which phase 
pairs are brought into contact. This turns out to be an im- 
portant avenue through which geometry can influence surface 
pattern. Furthermore, since the line tension between differ- 
ing phase pairs need not be the same, a three-phase system 
allows for an additional type of control parameter, one which 
can produce qualitatively different pattern behavior. 

In this work we analyze energy-minimizing structures of 
closed three-phase membranes. We observe that system fa- 
vors interface pairs which maximize differences in sponta- 
neous curvature. A pronounced effect of this preference is 
an effective attraction between phases of most disparate spon- 
taneous curvature, an indirect interaction driven by mutual re- 
sponse to geometric constraints. 

Fig- 03 illustrates the basic mechanism. Blocks of type /, 
77, III (red, green, blue) represent constituents with dif- 
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ferent preferred curvatures Cj, C'n, Cm. Possible four- 
interface arrangements are shown. The arrangement which 
maximizes the curvature difference at the boundaries is able 
to most closely accomodate its constituents' preferences when 
closed. 




Figure 1: (Color online) A cartoon representation of curvature-driven 
interfacial preference. Blocks of type /, II, III (red, green, blue) 
represent constituents with different preferred curvatures Cj, C'n, 
Cm. In this example \Ci - Cm\ > \Ci - Cn\ > \Cn - Cm\. 
Possible 4-interface arrangements of these constituents are shown. 
The arrangement which maximizes the curvature difference at the 
boundaries (bottom middle) is able to most closely accommodate its 
constituents' preferences when closed. 

We present cases of equal line tension and cases when line 
tension is pairwise preferential, showing that the geometry- 
induced attraction is robust enough to persist even in cases 
when it is most penalized by line tension. This study is fo- 
cused on characterizing the effective attraction. A full descrip- 
tion of the wide pallet of observed patterns will be published 
elsewhere 11911 . 



H. MODEL 

A. Continuum model 

We model the membrane as a two-dimensional surface. 
This approximation is appropriate for structures whose lat- 
eral dimensions are much larger than the membrane thickness. 
The system is considered in the strong-segregation regime, 
with each surface element identified with one of three phases. 
Each phase has a specific, predetermined spontaneous curva- 
ture, and interfaces between different phases are penalized 
with a line tension. The membrane is assumed to retain a 
closed spherical topology. We assume that the volume en- 
closed by the membrane is conserved, as is the area occupied 
by each phase. 



Our simulation searches for minimum energy configura- 
tions - that is, shape and phase arrangements which satisfy 
the area and volume constraints and minimize the total en- 
ergy. This energy is written as the sum of a bending term Fb 
and a line tension term Tg. 

The bending energy takes the form of a Helfrich functional 
ll2"oll . written as 

T b = J2j^( H -Cj) 2 dS 3 (1) 

j 

where j = 1,2,3 counts phases, H is the mean curvature, Cj 
is the spontaneous curvature of phase j, and k is the bend- 
ing rigidity. This term penalizes surface shapes whose local 
curvature deviates from its preferred value Cj . Since Cj de- 
pends on the local composition j, the bending energy couples 
shape and composition. Note that we do not give an index to 
bending rigidity k since we are assuming that the phases are 
equally bendable. 

Note also that the Helfrich functional typically includes a 
Gaussian curvature term weighted by saddle-splay modulus 
k. However, as with k, we are assuming that all k's are equal. 
Then, since the topology of the surface remains fixed, by the 
Gauss-Bonnet theorem 112 ill the Gaussian curvature term in- 
tegrates to a constant independent of the particular shape and 
therefore does not affect the energy. If the k's were not equal, 
the term would need to be retained. Phase-specific differences 
in bending rigidity k or saddle-splay modulus R p rovide addi- 
tional mechanisms for pattern formation [ 12, 22] and domain 
budding iH. 

In monolayers, the spontaneous curvature has an intuitive 
physical origin related to the properties of the head and tail 
groups [4]. Roughly speaking, molecules with larger effective 
head sizes and smaller effective tail sizes tend to pack into sur- 
faces of negative mean curvature, much as cones would pack. 
While the actual curvature properties depend on more than 
steric considerations, modeling lipid types with spontaneous 
curvatures has been shown to agree well with experiments 
112411 . In lipid bilayers, spontaneous curvature can arise from 
the effects of asymmetry between the inner and outer layers 
j25l - l27ll . Our model is therefore directly applicable to mono- 
layers, such as emulsions in which oil droplets are surrounded 
by multiple surfactant types, and to bilayers whose phases 
with distinct spontaneous curvatures are conserved, such as 
vesicles whose inner layer is uniform but whose outer layer 
is multiphase 112811 . To extend the model to general multicom- 
ponent bilayers would require treatment of separate bilayer 
properties. 

It is worth contrasting our deformation energy with that 
of solid-like membranes. There energy penalties associated 
with stretching and shear deformations can lead to buckling 
1I29I l3(ill and nonlinear conformation fluctuations ll3Tll in the 
homogeneous case, and in the two-component case elongated 
interconnected domains can arise [22, 32, 33]. In lipid vesicles 
with one fluid and one solid-like phase, complexes of stripes 
and polygonal domains have been observed experimentally 
113411 . These are outside the scope of our current study since 
none of our three phases include resistance to shear; it would 
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be interesting to explore interfacial preferences in these sys- 
tems. 

Line tension energy is given by ll35ll 



.(!) 



m, 



X tj d£ 
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where i,j = 1,2,3, \j is the line tension between phases i 
and j, d£ is a line element along i-j interfaces, and the line 
integral is calculated along the boundary 8y between phases 
i and j. This term favors segregation since it penalizes inter- 
faces between domains. 

Recently a novel class of molecules called linactants has 
been synthesized, which behave as 2-D analogues of surfac- 
tants: they migrate laterally in a fluid layer to boundaries be- 
tween immiscible phases, mediating between the immiscible 
constituents. Small amounts of linactant added to fluid layer 
have been found to reduce the line tension between immisci- 
ble phases ll3ql . This suggests that the line tension between 
phases can be viewed as a tunable design parameter. 

The surface area and volume constraints can be written as 
/ dSj = Aj (j = 1, 2, 3) and / dV^ = V , where Aj is the 
prescribed surface area of phase j, and Vq is the prescribed 
system volume. Note that there is an area constraint for each 
phase; the the system's total surface area and its composition 
fractions are both conserved. 

It is convenient to work with dimensionless parameters. 
Therefore, we measure energy in units of bending rigidity n, 
lengths in units of Rq, the radius of a reference sphere whose 
volume is Vq, curvatures in units of 1/Rq, and line tension in 
units of k/Rq. In these units, there remain nine independent 
input parameters: three spontaneous curvatures, three areas, 
and three line tensions. While this may seem like a daunting 
parameter space, we note that there is a great deal of symmetry 
and physical analogy, and we believe our choice of parameter 
pairs for simulation is representative and can provide useful 
insight into the phenomena at hand. 



B. Discrete model 

Our numerical method is based on a triangulation of the 
surface. Each vertex is identified with a phase type and sur- 
face area elements are represented by polygonal patches cen- 
tered at vertices. The discrete analogue of the bending energy, 
Eq. ((TJ, is computed as a sum over vertices, with the discrete 
mean curvature near each vertex v computed as ll37ll . 



^discrete ( v ) — 



i 

E \ A j 



(3) 



where i indexes edges for which v is an end point, |ej| is the 
length of the i th edge, tpi is the dihedral angle between two 
triangles sharing the i th edge, j indexes triangles containing 
v, and Aj is the area of the j th triangle. 

To compute discrete line tension, Eq. ©, we sum over 
all adjacent vertex pairs of different phase types rt35ll : 



E ^ij 

pairs of adjacent vertices 
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where denotes 



and j that have different phase 



types, cfj 1 and are the centers of the two triangles that 
have i and j as vertices, and my is the midpoint of their shared 
edge. Note that this differs from the approach used by Hu, et 
al. Ill 211 . as we take into account local triangle distortions when 
computing interface length. 

The surface area constraint is enforced by penalizing dif- 
ferences between A(i), the surface area associated with 
a vertex i, and the prescribed per-vertex value, F. d — 
§ — Ao(i)) 2 . The volume constraint is also enforced 

i 

by a harmonic potential penalty. 

Throughout all simulations, we set a = 2 x 10 5 and 
v = 10 4 . These values were tuned to ensure that, by the end 
of a run, the total volume and surface areas differed from their 
target values by less than 1%. Using <f>i to denote the compo- 
sition fraction of phase i, the prescribed surface area of phase 
i is 1.05 x 47r</>i. This is 5% larger than the area of a sphere of 
identical volume, excess being necessary to avoid the trivial 
case of an undeformable sphere. 

To find low-energy configurations, simulated annealing 
Monte Carlo simulations were performed following a linear 
cooling protocol. Three types of Monte Carlo moves were 
used: (i) a surface move, which attempts to perturb the posi- 
tion of a vertex by |Arj = 2.5 xlCP 3 , (ii) a phase-swap move, 
which attempts to exchange the phase type for two randomly 
selected vertices, and (iii) an edge flip, which cuts an edge 
shared by two triangles and reattaches it so that it spans the 
opposite previously-unattached vertices ll38l l39ll . All moves 
were accepted according to the Metropolis rules. 4 x 10 6 
sweeps, with a sweep defined as an attempted move of each 
vertex, were performed for each parameter set. Phase swap 
moves were performed every five and edge flip moves every 
ten sweeps. Initial configurations were random triangulations 
of a sphere, constructed from regular, Caspar-Klug triangula- 
tions ikOll such that vertex moves were constraint to a sphere 
and edge flips were performed at a very high temperature to 
ensure high acceptance rates. Any moves or edge flips that 
would result in unphysical edge crossings were rejected. The 
resulting configuration had a large number of vertices with co- 
ordination different than six. Finally, vertex types were ran- 
domly permuted. Each vertex was assumed to have a hard 
core of diameter l m in — 0.093, and each edge was endowed 
with a tethering potential with maximum length l max = 0.157 
such that I m ax /lmin = 1-688. These values were chosen to be 
tight enough to prevent membrane self-intersection but slack 
enough to allow edge-flips ll35[|4lll . 

Since the system's energy landscape is complicated, the 
computed configurations are not guaranteed to represent 
global minima. However, because independent runs start- 
ing from different random initial configurations reproduce the 
same qualitative features, we can regard the results as typical. 

Patterns are classified using graphs: each colored domain is 
identified with a colored node, and two nodes are linked if they 
share an interface. Then two configurations are considered 
to represent the same pattern if their graphs are isomorphic, 



respecting color. 



III. RESULTS AND DISCUSSION 



A sample of the varied patterns is provided in Fig. [2] where 
we present a phase diagram for a ^_ : _L : JL mixture of 
components I : II : III (red:green:blue). In this diagram, 
we vary Cm (spontaneous curvature of the blue phase) and 
A/,j7i (line tension between red and blue phases). We fix 
spontaneous curvatures Cj = 0, Cu = —1.0 and line ten- 
sions A/ // = Xjiju = 1.0. Results for each parameter set 
are classified according to pattern. We define pattern in terms 
of the number and arrangement of colored domains, irrespec- 
tive of shape. 
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Spontaneous curvature of blue (III) phase 
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Figure 2: (Color online) Phase diagram for composition ratio 7 : 
77 : 77/ — To : To : T5' varying spontaneous curvature Cm and 
line tension Xijii. All other parameters are fixed, with Ci = 0.0, 
Cu = —1.0, and \r,n = = 1.0. Phases I, II, III are 

shown in red (medium), green (lightest), and blue (darkest), respec- 
tively. Each example vesicle is shown twice: on the upper right, all 
phases are visible; on the lower left, one phase has been removed 
so that the pattern can be more clearly seen. Gray squares denote 
assorted structures with multiple green and red domains. All snap- 
shots were generated with the Visual Molecular Dynamics (VMD) 
package [ 42] and rendered with the Tachyon ray-tracer [43]. 

For present purposes, we note two features. The first is that 
a preferential line tension opens the possibility of qualitatively 
different pattern behavior and thus can be an important control 
parameter in the system's pattern selection. Furthermore, we 



observe a tranisition near Cm — —1.0, which is the dividing 
line where the spontaneous curvature of phase 777 switches 
from middle to largest in magnitude. 
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-2.8 -2.4 -2.0 -1.6 -1.2 -0.8 -0.4 
Spontaneous curvature of blue (III) phase 

Figure 3: (Color online) Analysis of a row of the diagram shown in 
Fig. [2] for A = 1.75. Composition fraction is set to : : 
and Cm is varied; Ci = 0.0, Cu = —1.0, Xiji = Xiijii = 1.0 
are kept fixed. Upper: Configurations with colored squares indicat- 
ing pattern grouping. Phases 7, 77, 777 are shown in red (medium), 
green (lightest), and blue (darkest), respectively. Middle: Length 
of each interface type to total interface length ratio r\ as a function 
of Cm. Interface fraction is represented by height between curves. 
Bottom: Spontaneous curvatures, with red, green, and blue lines cor- 
responding to spontaneous curvatures of phases of type 7, 77, and 
777, respectively. 

In Fig. [3] we analyze in detail the Xijii — 1.75 row from 
the diagram in Fig. [2] Note that I— III interfaces predominate 
in the region where Cm < Cu, even though such interfaces 
are nearly twice as costly as 7 — 77 and 77 — 777 interfaces. 

Fig.|4]examines a | : | : | composition ratio where all line 
tensions are set to 1.0, with (7/ = —1.6 and Cu = —1.0. In 
this row we observe that the spontaneous curvature of phase 
777 (blue) transitions from smallest in magnitude to middle to 
largest. Green-blue (77 — 777) interfaces are most favored on 
the left, where blue is largest in magnitude and green is small- 
est. In the middle, green (77) is smallest in magnitude and 
red (7) and blue (777) are close; here red-green (7 — 77) and 
blue-green (77 — 777) interfaces are about equally favorable, 
and red-blue (7 — 777) interfaces are disfavored. On the right, 
blue is largest in magnitude and red is smallest, and red-blue 
interface is most favored. Red and green are closest in magni- 
tude, and red-green interfaces are disfavored. The interfaces 
between phases of most disparate spontaneous curvature are 
most favored. 

In Figs. [5] and [6j we show that these interfacial preferences 
hold over wide slices of parameter space. Fig. |5]shows 7—777 
interface length as a fraction of total interface for the config- 
urations obtained from the parameter sets used in Fig. |2] In 
the regime where Cm < Cu < Ci, I — III interfaces pre- 
dominate until they become twice as energetically expensive 
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Red-blue (I-III) interface as a fraction of total interface 
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Spontaneous curvature of blue (III) phase 

Figure 4: Upper: Configurations with colored squares indicating pat- 
tern grouping. Phases I, II, III are shown in red (medium), green 
(lightest), and blue (darkest), respectively. Middle: Length of each 
interface type to total interface length ratio r\ as a function of Ci. 
Interface fraction is represented by height between curves. Lower: 
Spontaneous curvatures, with red, green, and blue lines correspond- 
ing to spontaneous curvature of phases /, //, and respectively. 
The composition fraction is set to | : i : i and Ci = —1.6, 
Cu — —1.0, — Xiijii — Xijii = 1.0 are kept fixed. 



Red-blue (I-III) interface as a percentage of total interface 




Figure 5: (Color online) Length of I — III (red-blue) interface as 
a fraction of total interface, for composition ratio I : II : III = 



j_ 

10 



j^j : jfi, varying spontaneous curvature Cm and line tension 



Xijit. All other parameters are fixed, with Ci = 0.0, Cu 
and A/./j = Xuju = 1.0. 



-1.0, 



as I -II or 77 - III. 

Fig. [6] plots I — 727 interface fraction in the system with 
| : | : | composition ratio, varying Ci and Cm. Xj jj = 
Xijh = Xjjjjj = 1.0 and Cu = —1.0 remain fixed. As 
expected, symmetry about the line Ci = Cm is evident. We 




Figure 6: (Color online) Length of / — /// (red-blue) interface as 
a fraction of total interface, for composition ratio I : II : III — 
| : j : |, varying spontaneous curvatures Cj and Cm. All other 
parameters are fixed. Cu = —1.0, and Xiji = Xr,m = Xnjn = 
1.0. 



observe the fraction of the / — III interface predominating 
in regions where Ci > Cu > Cm or Cj < Cu < Cm. 
Conversely, I — III interfaces are strongly disfavored in re- 
gions where \Cj - Cui\ < |C/ - Cu\ and \Ci - Cm\ < 
\Cu — Cm\. 

We point out that for large values of line tension or spon- 
taneous curvatures, well beyond the range used in this study, 
one observes interesting budding effects. Budding plays an 
important role in biological systems [44-46] and has been in- 



3SliZr|5Ql]. Effects 



vestigated in a number of studies [15, 
of interfacial preference on budding in a three-phase liquid 
membrane will be addressed elsewhere. 

In conclusion, we have seen that mutual response to ge- 
ometry acts as an effective attraction between phases of most 
disparate spontaneous curvature in a three-phase liquid mem- 
brane. This effect arises indirectly through the coupling 
of deformation and compositional arrangement, rather than 
through a direct component-component interaction. Nonethe- 
less, it is robust enough to compete with a countervailing line 
tension. In some cases it results in predominance of interfaces 
between phases least miscible by pure line tension considera- 
tions. Therefore, this system provides an example where ge- 
ometric constraints, rather than direct interactions, can dom- 
inate its conformation. Our findings suggest that an intricate 
interplay between the geometry and composition can lead to 
a rich phase behavior of complex fluid membranes. We hope 
that our results will stimulate further experimental and theo- 
retical work on these rich systems. 
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